Phase structures of strong coupling lattice QCD with overlap fermions at finite 

temperature and chemical potential 
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We perform the first study of lattice QCD with overlap fermions at finite temperature T and 
chemical potential \i. We start from the Taylor expanded overlap fermion action, and derive in the 
strong coupling limit the effective free energy by mean field approximation. On the (/x, T) plane and 
in the chiral limit, there is a tricritical point, separating the second order chiral phase transition 
line at small /x and large T, and first order chiral phase transition line at large and small T. 
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I. INTRODUCTION 

To study the nature of matter under extreme condi- 
tions in QCD is one of the most challenging issues in 
particle physics. Several novel phases have been sug- 
gested, such as quark-gluon plasma (QCP) 1] and color 
superconductivity 0- Precise determination of the QCD 
phase diagram on the (fJ>,T) plane will provide valuable 
information for these novel phases. 

Numerical simulation of lattice gauge theory (LGT) 
is the most reliable nonpcrturbative method based on 
the first principles. This approach has been successfully 
applied to the findings of the chiral and deconfmcmcnt 
phase transitions Q at finite T with zero [i. However, 
LGT experiences serious problems, like species doubling 
with naive fermions and complex action at finite fi. 

The Hamiltonian formulation of LGT at finite density 
El IE does not have the complex action problem. 
The complex action problem in Lagrangian formulation 
forbids numerical simulation at real fi. The recent years 
have seen enormous efforts @ 0, 0, on solving 
the complex action problem, and some very interesting 
information on the phase diagram for QCD at large T and 
small /i has been obtained. QCD at large /i is of partic- 
ular importance for neutron star or quark star physics. 

There have been several popular approaches to solving 
the species doubling problem of naive fermions. The stag- 
gered fermion approach preserves the remnant of chiral 
symmetry, but it breaks the flavor symmetry and doesn't 
completely solve the species doubling problem. The Wil- 
son fermion approach avoids the doublers and preserve 
the flavor symmetry, but it explicitly breaks the chiral 
symmetry; In order to define the chiral limit, one has to 
do nonperturbative fine-tuning of the bare fermion mass, 
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which seems to be an unnatural method from the physi- 
cal point of view|13|. 

These years have seen increasing interest in the over- 
lap fermion approach ^l^^], which is claimed to have the 
properties that chiral symmetry is preserved and species 
doubling problem may be solved 16J. Chiral symmetry 
in this approach is not the original form, and therefore 
one of the conditions of the Nielsen-Ninomiya (NO-GO) 
theorem^l is n °t satisfied. However, the Dirac operator 
of overlap fermion is nonlocal, and the computational 
costs [Tc^j for simulating dynamical overlap fermions are 
typically two orders of magnitude heavier than for the 
Wilson or KS formulations. It is also very tough to in- 
troduce the chemical potential into the action. Before the 
breakthrough of numerical algorithms for applying over- 
lap fermions to QCD thermodynamics, it is very useful 
to do an analytical study. At [i = and T = 0, this was 
made possible in Ref. , where the generalized overlap 
fermion approach was introduced by Taylor expansion of 
the Dirac operator. 

In this paper, we study the phase diagram of the strong 
coupling lattice QCD on the (/i, T) plane using Taylor ex- 
panded overlap fermions. In the chiral limit, we elucidate 
the phase structure and find lines of first and second or- 
der chiral phase transitions, as well as the tricritical point 
where the two lines join. 

The rest of paper is organized as follows. In Sec. 
ITT1 we shall give a brief review of the Taylor expanded 
overlap fermion approach, and rewrite the fermion ac- 
tion in terms of composite operators which transform co- 
variantly under the extended chiral transformation. In 
Sec. IIIII we introduce temperature and chemical poten- 
tial into the system and derive the effective free energy 
under mean field approximation. In Sec. IIVI we analyze 
the QCD phase diagram on the (fx, T) plane in the chiral 
limit. In Sec. the results are summarized. 
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II. FORMULATION OF GO FERMIONS 



It is easy to show that 



The action of lattice QCD is given by S = Sg + Sf,™, 
where 

Sg = -f X>(C/p + 02-2), 



Sf,w, = Sf + my^ j ip(x)ip(a 



(1) 



Here (5 = 6/g 2 , and U p is the ordered product of link 
variables U around an elementary plaquette. Sp is the 
massless fermion action: 



S F = a 4 ^2ip(x)D(x,y)ip(y), 

x,y 

with the overlap Dirac operator D(x,y) defined as^ r 

1 



aD = 1 + X- 



vxtx 



(2) 



(3) 



This operator is nonlocal and it is extremely difficult to 
do analytical calculations. An alternative is proposed in 
Refs. Pil|2(J, where the operator X is expressed as 



(4) 



with 



X = A- 










j 
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= S x+j,yUj(x) 
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(5) 



Then, one may expand D(x, y) in powers of the parame- 
ter t as[H H| 

aD x , y = 2d(A)5 xy + £.\2^ T J( x >y) + °( i2 )> ( 6 ) 

where we keep the contributions up to 0(t 2 ). In later 
discussions, we consider the case^ of negative A, which 
is expected to have the desired properties of QCD. 

In Ref. ^^j_it was shown that the Ginsparg- Wilson 
(GW) relation^ 



75D + D75 = aDj 5 D 



(7) 



implies an exact symmetry on the lattice, since the 
Nielsen-Ninomiya theorem is not applicable. It has been 
verified 0,113] that the i-expanded D(x,y) given by Eq. 
© satisfies the GW relation. 

The action is invariant under infinitesimal chiral trans- 
formation ip — ► ip + Sip and ip — > ip + Sip , where |16| 



Sip — £75(1 — aD)ip, Sip = £"075 



(8) 



(lp~75lp) -» (hsip) + M$ 1 - o o£) W>) 



(9) 



under the transformation JSj. This means that for m = 0, 
the vacuum has an exact chiral symmetry if the vacuum 
expectation value 



(r(l-ioD) r) =0. 



(10) 



Otherwise, the chiral symmetry is spontaneously broken 
if 



(11) 



These properties are very important for studying the chi- 
ral phase transition at finite temperature and chemical 
potential. 

Therefore, we could choose the quantity (ip(l — 
aD/2)ip) as the chiral order parameter, which will reduce 
to the conventional one in the continuum limit a — ► 0. It 
is more convenient to use the fermion fields q and q. They 
are related to ip and ip by [^3 



q = f. q = [ 1- j r. 



(12) 



The extended chiral transformation is defined by 

Sq = £759, sq = £<rr 5 . (13) 

The chiral order parameter is then given by 



Iqqr \l-\aD) 



(14) 



For convenience, we set the lattice spacing a = 1. Sub- 
stituting Eqs. 10, C21, and (JUJ) into Eq. (JTJ, we could 
rewrite the fermion action in terms of the new fermion 
fields q and q: 



Sf,™ = (l + —j^2q{x)D(x,y)q(y) 

x,y 

+ mY,<l(x)q(x)+0{t 2 ). (15) 



III. EFFECTIVE FREE ENERGY OF GO 
FERMION 

In the strong coupling limit, the gluonic action Sg 
vanishes, and S — > Sf,™,- We introduce the chem- 
ical potential \i into the action Eq. I|15f) by replac- 
ing the link variables in the temporal direction|22j with 

U 4 (x) -► e^Ui(x), and l/|(a;) -► e-^U\{x): 



S f 



1 + T ) ill ( £ u * ( x)q{ - x + ^ 

I 1 ^ a: j = l 



3 



- q(x + 3)ljU]{x)q(x)] 
+ ^[e"g»7 4 [/ 4 (i)g(i + 4) 

X 

- e-^q{x + 2)~f 4 Ul(x)q{x) 
+ m^2q(x)q(x). (16) 



Here d = 3 is the spatial dimensions. 

We use a notation x = (x, r) in which x and r rep- 
resents the spatial and temporal coordinate respectively. 
The temperature is given by T = N^ 1 with N T the num- 
ber of temporal lattice sites. We also use a representation 
of 7 matrices described in Ref. [23j . 

The partition function of the system is 

Z = J D[U 3 ]D[U A ]D[q]D[q]eK V {-SF, m ). (17) 

To derive the effective free energy F e fj = —hxZ/ 
we first integrate over spatial gauge link variables Uj, 
using the SU(N C ) group integration formulas and Taylor 
expansion: 



'«^(-( 1+ !)h£? 

x j — 1 

x [q(x)'y j U j (x)q(x +j) - q(x + J)j j uUx)q(x)' 



1-1 



2 fcY i 



A N c 



x ^2^2q(x)j jq (x + j)q(x + j)j jq (x), (18) 

X j = \ 

which should be a good approximation for Nf/N c < 3. 

Next, we linearize the four-fermion term in Eq. i|18|) 
by mean field approximation 6] 

d 

q(x)^q{x + j)q(x + jhjq(x) 

X j = l 

~ Y / d[Nf(v) 2 ~2vq(x)q(x)} 7 (19) 

x 

where v stands for the chiral condensate (qq). Substi- 
tuting Eqs. JJBJ and EH mto Eq- CZJ, the partition 
function becomes 

Z = [ D[U A ]D[q]D[q] exp \-S e " (U 4 , q, q; 0)1 (20) 



with 



S f,L (Ui,q,q;v) 

( 1 + y) ]§f E^^^+i) 



e^q(x + 4) 74 C/ 4 t (7(x)] 



J + — I — X!^/^ -2vq(x)q(c 



2 J N C A 2 
,^2q(x)q{x). 



(21) 



The hopping terms between q and q exist only in the 
temporal direction. 

To compute the remaining integrations, we make a 
Fourier transformation of the fermion fields 



1 

q(x,r) = -j= ^exp(ifc n r)g(f,n), 



n=l 



q(x, t) 



— Y exp(~ik n T)q(a 

' 1 V T 



(22) 



with kn = 2n(n — 1/2) /N T , and adopt the Polyakov 
gauge \ 



Ui(x,T) — diag 

= o. 

a=l 



e i0l(S) e i0 2 (S) e iSN c (x) 



(23) 



I.e., U^(x, t) is diagonal and independent of r. 

Substituting Eqs. and into Eq. J3UJ) and 

integrating out the Grassmann variables|25j. we have 

J D[q]D[q] exp £ ^Ag^ 



= (detA)^ 
where the determinant is 

N r N 



(24) 



det A 



nnn 

x a—1 n—1 

AC 2 / m\2 . 
"A 2 " I + 2" J (Sm n) 



, m\ d ( C 



1/2 



N, r 



nn 

x a — 1 



2C 



(1 + m/2) 



.14 

2cos[JV T (e o (5)-*^)] 



2 cosh[N T E] 



(25) 



with 



. 2(l + m/2)£(£) 2 0- 



E = arcsinhj 



mi + m/2) 



k n = k n + 9 a (x) - ifi. (26) 
We then perform the integration over U<± and obtain 

D[U 4 ](dctA) N f/ N - 



4 



n 



2C 
~A 



(l + m/2) 



4N f 



x 2cosh(iV c /iiV T ) + 



sinh [(N c + l)EN T ] 



Nj/N T 



(27) 



As a result, the partition function l|17l) becomes 

* = cx P ^-5:(i+w2) 2 ] g 5 ^(-) 2 ) 

_ _ /l AT 

2C, 



n 



.4 



(l + m/2) 



AN, 



2 cosh[iV c ^iV r ] + 



tix[{N c +l)EN T ]\ Nt/NT 



smh[EN T ] 



(28) 



Consequently, we obtain the effective free energy 
F eff = -lnZ/(£) 



m\ 2 dC 2 r T . .on 
1 + -) T^[N f (v) 2 ] 



2 J N r A 2 



NfT x In I 2 cosh 



- 42V/ In 



2C / m 



.4 




This formula is simple enough to make an analytical 
study of the vacuum, helpful in understanding the phase 
structure of strong coupling lattice QCD. 



IV. PHASE STRUCTURE IN THE CHIRAL 
LIMIT 

A. Chiral symmetry spontaneously breaking at 

H = T = 

In case of T = 0, the effective free energy (|2T))) reduces 
to a simpler form: 

T™ Feff = lv^' )2 " NcNf maX(/ "' E) ' (30) 



At /i = 0, the effective free energy becomes 

F eff [T = 0,a* - 0] = ^{v'f - N c NfE. (33) 

The rescaled chiral condendate v is determined by the 
conditions of minimizing the effective free energy 



dv' 

from which we get 



0. 



d{v) 2 



>o, 



' . \ 1/2 



(34) 



(35) 



with the spatial dimension and number of colors taken 
to be d = 3 and N c = 3 respectively. So the chiral 
condensate defined in Eq. (|14ll is 

v = (qq) = $ (l - ) = 0(t°)\ £ | ^ 0. (36) 

Therefore in the strong coupling limit, spontaneous sym- 
metry breaking of the extended chiral symmetry occurs 
at zero temperature and zero chemical potential. This 
qualitatively a gree s with that obtained from the external 
source method 19j. 



B. Chiral symmetry restoration at fi ^ and T = 

For small /x, the effective free energy is the same as Eq. 
i.e. 

F eff = N f v' 2 - W f E. (37) 

According to Eq. 134|) . the effective free energy reaches 
its global minimum F m i n = — 1.971.2V/ at v = 1.25. 

For large enough /i, the effective free energy H30|) 
reaches its minimum 



F min = -3JV//X 



(38) 



when v — 0, i.e., the system is in the chiral symmetric 
phase. 

Therefore, when \x increases from zero, the global min- 
imum changes from — 1.971^/ to —3Nf/i at the criti- 
cal chemical potential /i c , determined by — 1.971.2V/ = 
—3Nffi c . Therefore fi c = 0.657, at which the order pa- 
rameter v changes discontinuously, as shown in Fig. ^ 
This is a clear indication of first order chiral phase tran- 
sition. 



with 



and 



E = arcsinh|-^-u |, 



-' \ C -\ 
v = \^v\ 



(31) 



(32) 



C. Chiral symmetry restoration at fi = and 

In the small v region, the effective free energy takes 
the form 

F Bff (ji,T,v') = Cb(M,T) + Ca( M ,Ty 2 

+ C 2 ( f i,T)v^ + 0(v' 6 ), (39) 



0.0 - 

0.0 0.5 1.0 1.5 

Chemical Potential 

FIG. 1: Rescaled chiral condensate as a function of fi at T — 0. 
It drops discontinuously to zero at y, = fie = 0.657, which 
signals a first order phase transition. 
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FIG. 2: Rescaled chiral condensate as a function of T at /i = 0. 
It decreases continuously to zero at T = T c , which signals a 
second order phase transition. 



where 



C (fi,T) = -iV/Tln(2cosh(3M/T) + 4) 



Ci(jm,T) 



N f [ 1 



20 



5T + 2Tcosh(3/i/T) 



3T + 2cosh(3^/T) 

65-34cosh(3^/T) \ 
+ T 2 (5 + 2cosh(3M/T) J 



(40) 



According to the Landau's theorv|2q. a second order 
phase transition occurs when the coefficient of v' 2 be- 
comes zero. In the case fi = and T ^ 0, there is a 
critical value for T = T c at which Ci(0,T c ) = 0. In this 
case, v' vanishes continuously, implying a second order 
phase transition at the critical temperature T c , as shown 
in Fig. H 



D. Phase diagram on the (fi,T) plane 

Finally we study the phase structure on the (//, T) 
plane, based on above discussions. According to Eq. 
(|40|l , the coefficient of v' 2 vanishes when 



(rr s T c { 20-5T c 
/W ) = yarccosh I — — — 



(41) 



which corresponds to the second order chiral phase tran- 
sition line. When the coefficient of v' A becomes negative, 
first order chiral phase transition occurs at some positive 
value of Ci . 

Consequently, we have a tricritical point 



when C\ and Ci vanish simultaneously. These results 
indicate a second order chiral phase transition for T > 
Ttri and a first order one for T < Tt r i, as shown in Fig. 
The existence of the tricritical point is consistent with 
the recent result from the Hamiltonian approach with 
Wilson fermionQ- 




0.4 0.6 0.8 

Chemical Potential 



1.241, 



0.745. 



(42) 



FIG. 3: Phase diagram on the (/i, T) plane. The dotted and 
solid lines stand respectively for the second and first order 
transitions. The circle is the tricritical point. 



V. DISCUSSION AND OUTLOOK 

In the preceding sections, we have made the first at- 
tempt to investigate the phase diagram of Lagrangian 
lattice QCD with overlap fermions at finite temperature 
T and chemical potential /i. In the strong coupling limit, 
we discovered a tricritical point on the (/i, T) plane, sep- 
arating the second order chiral phase transition line at 
small fj, and large T, and first order chiral phase transi- 
tion line at large /j, and small T. 
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There are still some open questions. (1) The strong 
coupling is far from the continuum limit. (2) As in Refs. 

l2nl | , we considered only limited terms in Eq. (|19Jl • 
One may have to include infinitive terms in order to com- 
pletely suppress the doublers. Analytical calculations 
will be extremely difficult 0, |2£j], even at /i = T = 0. 
Further study of QCD with overlap fermions at finite 
temperature and chemical potential, ether analytical or 
numerical, will be very interesting. 
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